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ABSTRACT 

In hierarchical models of structure formation, the time derivative of the halo mass 
function may be thought of as the difference of two terms - a creation term, which 
describes the increase in the number of haloes of mass m from mergers of less mas- 
sive objects, and a destruction term, which describes the decrease in the number of 
m-haloes as these merge with other haloes, creating more massive haloes as a re- 
sult. The first part of this paper focuses on estimating the distribution of times when 
these creation events take place. In models where haloes form from a spherical col- 
lapse, this distribution can be estimated from the same formalism which is used to 
estimate halo abundances: the constant-barrier excursion-set approach. In the excur- 
sion set approach, moving, rather than constant-barriers, are necessary for estimating 
halo abundances when the collapse is triaxial. First we generalise the excursion-set 
estimate of the creation time distribution by incorporating ellipsoidal collapse. Then 
we show that these moving-barrier based predictions are in better agreement with 
measurements in numerical simulations than are the corresponding predictions of the 
spherical collapse model. In the second part of the paper we link the creation times 
distribution to the creation term mentioned above. For this quantity, the improve- 
ment provided by the ellipsoidal collapse model is more evident. These results should 
be useful for studies of merger-driven star-formation rates and AGN activity. We also 
present a similar study of creation of haloes conditioned on belonging to an object 
of a certain mass today, and reach similar conclusions - the moving barrier based es- 
timates are in substantially better agreement with the simulations. This part of the 
study may be useful for understanding the tendency for the oldest stars to exist in 
the most massive objects, and for star formation to only occur in lower mass objects 
at late times. 
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INTRODUCTION 



In hierarchical clustering models, self-bound dark matter 
haloes increase their ma s s by merging with other haloes 
llPress fc Schechterl 1 19741 : iBond et all Il99ll : lLacev fc Coll 
Il993l h These mergers are expected to a ffect the galaxy pop- 
ulations hosted by th e merging haloes (|White fc ReesHl978l : 
IWhite fc Frenklll99llh possibly trig gering star formation or 
AGN activity dEfstathiou fe Rees! Il988l ; lHaehnelt fc Reesl 
1 19931 : iKauffmann fc Haehneltll200d h As a result, there has 
been some interest, both analytic and numerical, in estimat- 
ing when such mergers happen. 

Let n(m\t) Am denote the number density of haloes with 
mass in the range dm about m at time t. As a result of 
mergers, An/ At is the sum of two competing effects - the 
number of objects in a given mass bin increases if smaller 
mass objects merge to form an object of precisely this mass 

* Email: jmoreno, shethrk@physics.upenn.edu carlo.giocoli@unipd.it 



- an event we call creation - or the number decreases as 
objects of this mass merge with others, thus depleting the 
number in the bin - an event we call destruction. Thus, the 
time derivative of the halo mass function is the difference of 
the creation and destruction rates: 



An/ 'At = C(m, t) - D{m, i). 



(1) 



For the star formation and AGN problems above, one is 
more interested in the creation rates C(m,t) than in An /At. 

Given n(m\t), it is easy enough to take the time deriva- 
tive; the problem is to separate An/ At into its two contribu- 
tions. Roughly speaking, low mass objects may have under- 
gone significant mergers in the past but they are not being 
created in merging events any more - their evolution is ex- 
pected to be dominated by the destruction term. In contrast, 
extremely massive objects are undergoing substantial merg- 
ing activity at the current time, and the time derivative of 
the halo mass function should be a good estimator of the 
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creation rate of these objects. But quantifying the general 

case requires a richer mo del. 

Early work used the lPress fc Schechterl (|l974T j form for 
n(m\t), and advocated equating the 'positive' term in dn/dt 
with the creation r ate, and the 'negative ' term with the de- 
struction rate (e.g., Ha ehnelt et al, . 1998). But this is clearly 
not a solution at all, since it provides no rule for how 
to determine what one should correctly equate with 'pos- 
itive'. For example, if dn/dt = V — M , there is no par- 
ticular reason why one could not have written the right 
hand side as (V — e) — (TV — e). The second part of this 
paper is de voted to extract i ng the cre ation r a te C from 
dn/dt. See iBlain fc Longairl ljl993al lbh. ISasakil (|l994h and 
iKitavama fc Sutol (|l996h for other attempts to solve this 
problem. 

Before addressing the creation term C (the number 
density of mergers per Gyr), in the first part of the pa- 
per we discuss the distribution of times c(t\m) when these 
creation events take place (i.e., c = C/ j C dt is the rate 
normalised by the total number of creation events that will 
ever occur). Although c and C differ only by normalisation 
constant, it turns out that c is somewhat easier to model. 
This is because the excursion set formalism fro m which the 
Press - Schechter mass func tion can be derived ijBond et alJ 
Il99ll : lLacev fc Coll Il993h . carries with it a prescription 
for computing c( t \m), the distribution of creation times 
IPercival fc Miller! l|l999h . In this case, the functional form 
of c(t\m) is very similar to that of f(m\t) = mn(m\t)/ p, 
where p is the mean comoving background density. 

Since that time, interest has shifted to functional forms 
for n(m\t) which more closely approximat e the abundances 



teractions include Silk fc White! (Il978l). ICavaliere et al 



Jl991aT), ICavaliere et all (Il991bl), ICavaliere et al] (|l99: 



?. IShetii fc Pitmanl ij 19971 ) andlMenci et al] 200^) For a 



meas u red in numeri c al simulations (e.g . . Shcth fc Tormenl 



1999; Warren et al. 2006; Rccd ct al. 2007; Lukicetal 



20071 : iTinker et alll2008n . So it is interesti ng to ask how 



the cr eation time distributions are modified. IPercival et al] 
(2000) argue that the relation between the functional forms 
of c(t\m) and f(m\t) should survive these modifications, and 
show that this does indeed provide a good description of halo 
creation in simulations. However, although they use intu- 
ition from the excursion set approach to motivate their argu- 
ments, their method side-steps the generalisation of the ex- 
cursion set approach from which the modified mass functions 
may be derive d - this is the ellipsoidal collapse 'moving' bar- 
rier a pproach (|Sheth. Mo fc Tormenll200ll ; ISheth fc Tormenl 
2002). The first goal of this paper is to calculate the creation 
time distribution self-consistently within the moving barrier 
excursion set approach. We do find that c and / are sim- 
ply related, but that this is actually extremely fortuitous - 
Appendix IA3I demonstrates that this scaling does not hold 
generally. 

Getting the normalisation constant which relates c to 
the creation rate C (m, t) is a more challenging problem. 
IPercival et ail l|2000l i obtained this quantity by matching 
the creation time distribution to the rate measured in 
N-body simulations in the low redshift regime (see also 
IPercival fc MilleU Il999l ). but they acknowledge that they 
have no theory for the normalisation factor. One pos- 
sible solution to this problem is to explore the evolu- 
tion of the halo population in terms of coagulation the- 
ory, where the cre ation and destructi o n ter ms are esti- 
mated separately (|Smoluchowskil Il916l . Il917l ). Early ap- 
plications to galaxy formation and dark-matter halo in- 



more recent treatment, see lBenson et al] "( 2005 ) and lBensonl 
|200Sft . 

For white-noise initial conditions, both the Smolu- 
chowski and the Press-Schechter excursion-set expressions 
for n(m\t) agree, so, f or this case, the creatio n and destruc- 
tion rates are known (jSheth fc Pitmanl fl997h . However, ob- 
taining the rates for more general initial conditions, or for 
the modified mass functions that are of more current inter- 
est, remains unsolved. The second goal of the present paper 
is to provide a model for the creation rate of dark matter 
haloes that is informed by both coagulation theory and the 
modified excursion set approach with moving barriers. 

We also study the problem of how halo creation is 
modified if it is known that the merging haloes are bound 
up in objects of mass M at some later time T. The ex- 
cursion set theory provides a way to compute the condi- 
tional mass function N(m\t, M, T). We show how the condi- 
tional distribution of creation times c(t\m, M, T) is related 
to f(m\t,M,T) = (m/M)N(m\t,M,T). For the conditional 
rate, the problem is to separa t e dN /dt into creation and 
destruction components. Sh ethl l|2003l ) argues that this con- 
ditional distribution may be the basis for und erstanding the 
phen omenon known as down-sizing (also see iNeistein et al] 
l200fj) . 

Section [2] provides a brief summary of the excursion set 
approach and shows why the creation times associated with 
different barriers differ. Sections [3] and U] compare the pre- 
dicted and measured N-body creation time distributions and 
creation rates, respectively. Conditional versions of these 
quantities are included in the appropriate sections. A final 
section summarises our findings. Technical details of the cal- 
culations are provided in the Appendices. 

A few final remarks regarding the different uses of the 
term 'halo creation' are in order. The first part of the paper 
focuses on the creation time distribution, c(t\m), which can 
be derived within the excursion-set formalism. The second 
part focuses on the creation rate, C(m,t), the first term in 
the coagulation equation (equation [1]). The former is a nor- 
malised time distribution, while the latter is not normalised. 
(The normalised distribution is denoted with lower case c, 
while the un-normalised rate is denoted with capital C ) An- 
other source of confusion is that halo 'creation ' is di stinct 
from halo 'formation'; following lLacev fc Cokl l|l993h . the 
latter is typically defined as the t ime that an object first 
reaches half its current mass. See iGiocoli et all (|2007h for 
an explicit calculation showing how creation and formation 
are related. 



2 EXCURSION SETS AND MASS HISTORY 

In the excursion set approach, the problem of estimating the 
halo abundances is mapped to one of estimating the distri- 
bution of the number of steps a Brownian-motion random 
walk m ust take before it first crosses a barrier of specified 
height (|Bond et al.lll99lh . In this approach, the height of 
the barrier plays a crucial role. The Press-Schechter mass 
function is associated with barriers of constant height - 
such barriers arise naturally in models in which haloes form 
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Figure 1. The mass history associated with a random walk (jagged line) depends on the properties of the barrier. Panels show constant 
(left) and square-root (right) barriers. Time increases as 8 decreases and mass decreases as S increases. The filled circles on the random 
walk denote the history (red for spherical and blue for ellipsoidal collapse), and the horizonal jumps denote mergers. In the right panel, 
the red circles were kept to emphasize that different barriers predict different mass histories for a given random walk. For reference, the 
horizontal dotted line (curve) denotes the barrier associated with the present. 



from a spherical collapse model. In constrast, the more ac- 
curate mass functions may be related to barriers whose 
height increases monotonically with the number of steps - 
such barriers arise naturally in ellipsoidal collapse models 
|Sheth. Mo fc TormenlkoOlf). 

Following ISheth fc Tormenl l|2002h . we will be interested 
in barriers of the form 



B(S,S C ) = y/qS c il + p 



S 



qSi 



(2) 



Here 5 C is the overdensity required for spherical collapse - 
it is a monotonically decreasing function of time, given by 
5co/D(t), where <5 c o — 1.686 and D(t) is the growth fac- 
tor. S is a monotonically decreasing function of halo mass, 
given by a 2 (m), the variance of the initial density fluctua- 
tion field. The Press- Schechter mass function is associated 
with (q,/3, 7) = (1,0,0), whereas ellipsoidal collapse has 
(0.707,0.47,0.615). 

When (3 = (the barrier associated with spherical col- 
lapse) all walks are guaranteed to cross the barrier, and the 
barriers associated with two different times do not intersect. 
However, if j3 > and 7 > 1/2, e.g., for the ellipsoidal col- 
lapse barrier, then not all walks cross the barrier, and the 
barriers associated wit h two different times may intersect. 
ISheth fc Tormenl {2002) suggest that the intersection of bar- 
riers may allow one to represent the possibility that haloes 
can fragment. This is problematic when discussing the cre- 
ation rate problem, which assumes that fragmentation never 
occurs. 

For this reason, we study the limiting case of 'square- 
root' barriers for which 7 = 1/2: 



This family of barriers is particularly interesting because an 
analy tic solution to the first crossing distribution is avail- 
able (Breiman 1966), although slightly cumbersome (see Ap- 
pendix [XT}. Because 7 = 1/2 is not very different from the 
value associated with ellipsoidal collapse, one might have 
expected the predicted halo abundances associated with 
square-root barriers to provide a reasonable description of 
simulations. We show below that this can be achieved if one 
sets (q, j3, 7) = (0.55,0.5,0.5) (see Figure [3}. Halo merger 
and formation histories associated with this model are also 
in good agreement with simulations (|Giocoli et alj 120071 ; 
iMoreno et. ai~ll2008T i. 

The dependence on 5* of the square root barrier means 
that it is more like the ellipsoidal than spherical collapse 
barrier (which has constant height, independent of S). How- 
ever, there is one important respect in which the square root 
model is very like the constant one. Consider the barriers as- 
sociated with two different times. For square-root barriers, 
the difference between the barrier heights is 



B(S,5 c2 )-B{S, 6*) 



(S c2 
S C 2 - 



+ 0VS) - (<5ci + pVs) 

-Sci. (4) 



B(S,S C ) = ^s c +pVs. 



(3) 



Notice that this difference is independent of S. This is also 
(trivially) true for constant barriers, but it is not true for any 
other values of 7. In this respect, the excursion set model 
based on square-root barriers is extremely special. This will 
be important later. 



2.1 Mass history for different barriers 

Figure [T] illustrates the relation between Brownian motion 
random walks and the mass growth history of an object. 
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Consider first the panel on the left. The jagged line shows 
an example of a random walk — this walk represents the 
run of smoothed overdensity around a randomly chosen po- 
sition in the initial fluctuation field, as the region over which 
the overdensity is smoothed changes from large (left) to 
small (right). (The plot actually shows the initial overdensity 
evolved to the present time using linear theory — it differs 
from the initial overdensity by a multiplicative constant.) 
The initial overdensities are all small compared to unity, so 
one may associate a mass with each smoothing scale: this 
mass is larger for the larger smoothing scales. 

Consider a horizontal line, and consider the places 
where it first intersects the random walk, as the height of 
this line, this barrier, is raised. Clearly, this position shifts 
to the right as the barrier is raised (red filled circles) — mass 
decreases as redshift increases. Whereas the halo abundance 
problem corresponds to fixing the barrier height 8 C (to illus- 
trate, the height of dotted line corresponds to 8 c o = 1.686) 
and asking for the distribution of S values at which the bar- 
rier is first crossed, the halo creation problem corresponds 
to asking for the distribution of S c values for a fixed S. We 
will use f(S\S c ) dS to denote the first crossing distribution, 
and c(S c \S) dS c to denote the distribution of creation times, 
where, for haloes of a given mass m, c(S c \S) dS c = c(t \m) dt. 

Notice that the mass increases relatively smoothly 
sometimes, and rather abruptly at others. For instance, in 
the interval 4.2 < S < 4.8 in the left panel, mass decreases 
smoothly. Compare this situation to the sudden jump from 
S ~ 4.8 — » S ~ 8.7 (red long-dashed line). Thus, this walk 
does not contribute to the creation time distribution for any 
values of S between 4.8 and 8.3. But it does contribute in 
the calculation of halo abundances for every <5 C . 

It is interesting to compare this mass accretion history 
with that shown in the right panel. The same walk is shown 
in both panels. However, now the horizontal dotted line at 
Sco has been replaced by a curve that has 5-intercept at 
^/qS c o, and increases with increasing S — this is the square- 
root barrier (equation [3) associated with the same epoch 
as the constant one. The comparison clearly shows that the 
mass accretion history (blue solid circles) depends on the 
barrier shape. For instance, the values with S < 4.3 and 
S ~ 8.8, 9.7 and S < 10.4 are no longer included in the 
mass accretion history (red solid circles). Moreover, even 
if a point on the random walk happens to be part of the 
mass history in both cases, its associated time is different 
under the two barrier prescriptions. To illustrate, the point 
at S ~ 4.4 is at the present for the square-root barrier case 
and in the! past when constant barriers are used. As a result, 
the halo mass function and the creation-times distribution 
are modified. One of the goals of this paper is to quantify 
these changes. 

2.2 Self-similarity 

Barriers of the form given in equation |2]) are self-similar, 
in the sense that if S c is increased by a factor ft, then so is 
y/S. As a result, the first crossing distribution f(S\S c ) and 
the distribution of halo creation times c^lS) are both sim- 
ply functions of S^/S. In the present paper we will denote 
this variable as v if S c is fixed and as v c if S is fixed. In 
other words, if equation ([2]) describes ellipsoidal collapse, 
then f(S\5 c )dS = f{u)du and c(5 c \S)dS c = c(v c )dv c . Ap- 
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Figure 2. The creation time distribution in self-similar form. The 
variable u c denotes at fixed mass. Triangles and circles show 

the distribution measured from an ensemble of random walks with 
constant and square-root barriers respectively. Dotted and solid 
(dashed) lines show the associated pred ictions. The long dashed 
curve shows the result of inserting the ISheth fc Tormenl (Il999h 
form for f(u)du into equation {5J. 

pendix |A"T1 shows that 

c(f c ) dv c = A^/i^f(v c ) dis c (5) 

for the constant (.4 = and square root (.4 ~ 2) bar- 

riers. This simple relationship is one of the central results of 
this paper, as is the warning that it does not hold in gen- 
eral. Appendix A uses a system of linear barriers to illustrate 
when this simple result does not apply. 

To test equation © for constant and square-root barri- 
ers, we generated 10 5 random walks with 10 4 steps between 
S = and S = S(m p ) ~ 28. Then we stored the correspond- 
ing mass histories for the constant and square-root barriers 
(e.g., solid circles in Figure [TJ. Every (5*, <5 C ) point along 
the history has a corresponding v c . To study how the cre- 
ation time distribution depends on S, we could have chosen 
the subset of walks which have the correct value of S, and 
then found the distribution of v c = values for those 

walks. However, the self-similar scaling above means that 
c(i/ c ) should be the same for all S. As a result, there is no 
need to select a subset in S before binning in u c . If we simply 
bin up all the v c values, whatever the associated values of 
S, then we can compare the result with the predicted c(u c ). 

The symbols in Figure [2] show the creation time dis- 
tributions for the constant (triangles) and square-root (cir- 
cles) barriers: the distribution associated with the square- 
root barrier is broader and peaks at slightly higher v c . The 
curves show the predicted creation time distributions (equa- 
tion 0; they are in excellent agreement with the measure- 
ments. Equation ([5]) shows that these creation time distri- 
butions depend on the functional form of the first cross- 
ing distribution f(v). For the case of square-root barriers, 
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Figure 3. Comparison of the GIF2 mass function (symbols) with that derived from the first crossing distribution of a square-root barrier: 
equation ((3} with (g,/3,7) = (0.55,0.5,0.5). The constant barrier Press-Schcch ter result is p lotted for completeness. Bottom panel shows 
the ratio of both data and theory curves to the functional form of lSheth &: To rmen (199^). 



we show the iBreimanl l|l966l ) exact but complicated expres- 
sion for f(v), and a muc h simpler approximation for it from 
ISheth fc Tormenl {2002). The two curves are almost indis- 
tinguishable. 

The use of barriers which scale self-similarly (equa- 
tion [2]) was motivated by the observation that, when ex- 
pressed as a function of v, halo ab undances in simulation s 
could be scaled to a universal form (|Sheth fc Tormenll 19991 ). 
We have added the long-das hed curve in the F i gure, which 
shows the result of using the lSheth fc Tormenl (jl999T ) func- 
tional form for f(v c ) in equation ([5]); it is almost indistin- 
guishable from the curves associated with the square-root 
barrier. 



surements, we now study if they provide a good description 
of halo creation in N-body simulations. 

We u se data, f rom the GIF2 sim- 
ulation l|Gao et al.l 120041 ') . available at 
http://www.mpa-garching.mpg.de/Virgo , which 
followed the evolution of 400^ particles of mass 
m p = 1.73 x lO 9 fr _1 M0 in a box of size of 
110/i _1 Mpc, in a flat ACDM cosmology with parame- 
ters (n m ,a 8 ,h,tt b h 2 ) = (0.3,0.9,0.7,0.0196). Haloes were 
identified at 50 outputs equally spa ced in log 1fl ( l + z ) 
between 1 + z = 20 and 1 + 2 = 1. See lGiocoli et all (|2008t) 
for more details about the post-processing of the simulation. 



3 CREATION TIME DISTRIBUTION 

In Section[2]we discussed the creation time distribution from 
the excursion set point of view. Having shown that the ana- 
lytic expressions accurately reproduce our Monte Carlo mea- 



Haloes were labelled as having been 'created' if at least 
half of their particles were not observed in a more mas- 
sive halo at an earlier time. This is essentially the method 
adopted by IPercival fc Milled ljl999l ) - and we refer the in- 
terested reader to that paper for details. 
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Figure 4. Distribution of creation redshifts for a number of bins in halo mass. Filled circles show measurements in the simulations, and 
open triangles and circles show analogous measurements made from sampling our constant and square-root barrier random walk ensembles 
similarly to the simulations. Dotted curve shows the prediction associated with a constant barrier; the exa ct square-root solution a nd its 
series approximation are the solid and dashed curves; long-dashed curve shows the result of inserting the lSheth fc Tormenl jl999h form 
into equation {5jl. 



3.1 Distribution of creation redshifts 

The creation time distributions we measure in simulations 
and shown as filled circles in Figure [4] are normalised to 
unity. However, the simulations only sample <5 C at epochs 
before the present time, whereas the theory curves assume 
that ^ S c ^ oo. Therefore, for haloes of mass m, we set 



c(z\m) dz 



c(S c \S„ 



f g °°d6> c c(5> c \Sr, 



d5c 
dz 



dz, (6 C > 5 c0 ). (6) 



These are the curves in Figure [4] The open triangles 
and circles show the creation times measured in constant 
and square-root random walk ensembles sampled at the 
same redshifts as the simulations. We studied bins of size 
dlog 10 m = 0.2 in mass centred at log 10 (m/m*) = 0.5 
to —3 in steps of —0.5. It is common practice to express 
halo masses in terms of the typical mass m*(z), defined 
by S(m i< (z)) = 5't(z). Throughout this work, m* (with 
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Figure 5. Universality of the distribution of halo creation times. The same theory curves as in Figure [2] are shown. Different symbols 
show results for different halo masses in the G1F2 simulation data, as indicated. 



no jz-dependence) denotes m*(z 
m + = 8.7 x 10 12 M o /i _1 ~ 5030m D 



0). In this cosmology, 



A couple of remarks are in order. Only the haloes with 
the highest redshift in each mass bin were treated as newly- 
created. These measurement were tested with different bin 
sizes (not shown), yielding similar results. One limitation 
is that if dlog 10 m is too small, most bins are empty, and 
the data does not follow a continuous curve. One should not 
take dlog 10 m to be too large - in particular, the mass bins 
for the different values of m should be disjoint. Our choice 
of dlog 10 m satisfied both criteria. Lastly, we found that our 



choice of redshift bin dlog 10 (l + z) = 0.05 was sufficiently 
large to capture enough creation events and sufficiently small 
for comparison with the different theory curves. 

The solid circles in the Figure denote the N-body mea- 
surement, where only haloes with m > 10m p are consid- 
ered (notice that the lowest-mass panel has no black filled 
circles). In all cases, improvement over the location of the 
peaks is seen when the square-root barrier is used. How- 
ever, these curves are slightly broader than those traced by 
the simulation data, making the height of these normalised 
distributions lower. 
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Figure 6. Conditional distribution of creation redshifts. Symbols and style as in Figure l4l We plot m = m*/10 (left panels) and 
m = m*/100 (right panels) conditioned to end up in haloes of mass M = M+ (top panels) and M = 10m* (bottom panels). In all cases, 
T denotes the present time. 



3.2 Self-similarity of halo-creation 

The excursion set model suggests that if the halo mass func- 
tion f(v)dv can be scaled to a self-similar form, then the 
creation time distribution c(v c )dv c is also self-similar. To 
test this we have scaled the values of S c (z) associated with 
each mass bin in the simulations to v c — S c (z) 2 /S(m), and 
measured the resulting distribution of u c . However, because 
all mass bins sample the same range in 8 C , they sample dif- 
ferent ranges in v c . We account for this by dividing the mea- 
sured distribution of u c by a normalisation factor given by 



JJ° du' c c(u' c ), where v c o = Sco/S(m) and c(v c ) is associated 
with the lSheth fc Tormenl l| 19991) formula. 

Figure [S] shows the result: different symbols show the 
rescaled distributions associated with the various masses. 
Note that they do indeed appear to trace out a universal 
curve. The various smooth curves show the constan t bar- 
rier, square-root barrier, and lSheth fc Tormenl (|l999l ) based 
predictions. The symbols approximately split the difference 
between the constant and square-root barrier models. 
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3.3 Conditional distribution of creation redshifts 

So far we have discussed the unconditional creation-time 
distribution c(t \m) and its relation to f(m\t). In this section 
we study the creation time distribution, c(t \m, T, M), of m- 
haloes at time t conditioned to be bound up in M-haloes at a 
later time T. We also discuss how this quantity is related to 
f(m\t, M, T), the fraction of mass in m-progenitors at time 
t of a final halo of mass M at time T. The latter is derived 
from the excursion-set theory by setting 

f(m\t,M,T)dm= f(S\5i,S ,S )dS, (7) 

where S = S(m), S = S(M), S 1 = S c (t) and S = S C (T). 
The right-hand term is the crossing distribution of a barrier 
B(S,8i) by random walks with origin at (So, B (So, So)). In 
Figure [l] this amounts to shifting the origin from (0, 0) to 
(So, So) (left panel) or to (So,^/q5o + fty/So) (right panel). 
In this work, the conditional mass function is denoted by 
N(m\t, M, T)dm = (M/m)f(m\t, M, T)dm. 

In essence, conditioning is equivalent to finding the (un- 
conditional) crossing distribution of a barrier 

B(s,S 1 ,So)=B(s + So,S 1 )-B(S ,So), s = S - So- (8) 

In the constant barrier problem, B — S\ — So- This means 
one simply replaces all the unconditional expressions given 
previously with S c (z) — > Si — So and S —> s = S — So- As a 
result, the only change occurs in the s elf-similar variable v 
(and u c ): S 2 JS -> (Si -So) 2 /(S - S o) jLacev fc Cole! 1 19931 ; 
IPercival fc Miller! 1 19991 : IShethll2003l ) . 

The square-root barrier is slightly more complicated: 

B(s, o c ) = a c + PVs + S , a c = Si- 8 - p^/S~o~. (9) 

Because equation © is not quite the same functional form 
as equation ©, the first crossing distribution is not simply 
a suitably rescaled version of the unconditional distribution. 
Rather, it is a function of 77,3 = a c /y/~So and s/So- Never- 
theless, the logic one follows to arrive at the creation time 
distribution is the same. In particular, the conditional ver- 
sions of c and / are simply related: 

c(r h3 \S/So)d V p = A {s/ s )f(S/So\vp)d(S/S ), (10) 

where now the A factor is a function of S/So- This is an- 
other central result of this paper (compare to equation [5]) . 
Expressions for the exa ct solution and the corresponding 
ISheth fc Tormenl (|2002h approximation are given in Ap- 
pendix [A5] 

Figure [6] shows the conditional distribution of creation 
redshifts for m/m* = (0.1,0.01) that end up in haloes of 
mass M/m* = (1, 10) today. Filled symbols show the GIF2 
measurements, open circles show our Monte Carlos with 
square root barrier, and open triangles show Monte Carlos 
with a constant barrier. Smooth curves show the correspond- 
ing predictions - notice that they are in excellent agreement 
with the Monte Carlos. 

The simulation bin sizes dlog 10 m and dlog 10 (l + z) were 
used as in the unconditional case. Selecting haloes to be 
conditioned to belong to a final M-halo reduces the number 
of creation events significantly. We selected haloes bound to 
end-up in clumps with mass in a bin of size dlog 10 M = 0.5. 
As in the unconditional case (Figure the distributions 
peak at higher redshifts and are slightly broader for lower m 
(compare left and right panels) . The same trends are seen as 
one increases the final M (compare top and bottom panels). 



In general, the moving barrier based curves provide a 
much better description of the simulations, although the 
agreement is by no means perfect. For example, the square 
root barrier tends to produce distributions which are slightly 
broader than those in the N-body simulation. This effect is 
more pronounced i n the right-hand side panels. A similar 
effect was found bv lMoreno et. al~l ()2008h for the formation 
time distribution, which is related to (but diffe rent from) the 
creat ion time distribution of interest here (see iGiocoli et al.l 
120071 . for details). In that paper, we speculate that the dis- 
crepancy there was due to non-M a rkovia n effects. We refer 
the interested reader to lPan et al.l l)2008h for a discussion of 
this topic. 



4 HALO CREATION RATES 

Extracting the creation rate from dn/dt (or dN/dt) is a 
non-trivial problem. In this section we make use of halo co- 
agulation theory to estimate this quantity. 



4.1 Unconditional rate 

In this formalism, the halo mass function n(m\t) obeys 

^l=C(m,t)-D(m,t), (11) 

where 

„. . f m K(m! , m — m';t) , ,. , , /, , , / 

C(m,t) — I — -n(m \t)n(m — m \t)dm , 

Jo 2 

(12) 

is our creation term, and the destruction term is 

D(m,t) — I K(m,m';t)n(m\t)n(m'\t)dm' (13) 
Jo 

(|Smoluchowskilll917l ). In these expressions, the coagulation 
kernel K(m,m';t) is symmetric in m and m' . 

Few analytic solutions to Smoluchowski's equation ex- 
ist. However, when the kernel is additive in mass, then the 
associated mass function is given by the Press-Schechter for- 
mula for white-noise i nitial conditions (|Silk fc White! IT978I ; 
ISheth fc Pitmanlll997l ). Of course, white-noise is a bad ap- 
proximation to the initial conditions in the CDM models 
of current interest. Moreover, we have shown that ellip- 
soidal collapse gives a better description of halo abundances 
and creation times than Press-Schechter (spherical collapse). 
Nevertheless, the expression obtained for the creation term 
in that special case, 

C(m,t) = pmn(m|t)|^p| (14) 

ijSheth fc Pitman! Il997h . will serve as a guide (see Ap- 
pendix B for more details regarding coagulation with spher- 
ical collapse and white-noise initial conditions). 

This prescription has two interesting properties. First, 
it is related to the creation time distribution in a simple 
way: 

C(m,t) = g(m)c(t\m) (15) 
(please see Appendix |B2}. This is equivalent to saying that 
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Figure 7. Halo creation rates. Symbols, line styles and choices of mass as in Figure [4] As in Figure [I] the lowest panel on the right 
shows not data because this case is beneath the resolution of the N-body simulation. 



one obtains c by normalising C: 

C(m,t) = g(m)c(t\m) _ g(m)c{t\m) = c n\ m ) 
fC(m,t)dt J g(m)c(t\m)dt g(m) J c(m,t)dt c{ - < m > 

(16) 

This indicates that all the time dependence in the creation 
rate is encoded in the creation time distribution, which was 
amply studied in the first part of the paper. In particular, 
in Section [3] we showed that the mass function and the cre- 
ation time distribution are related in a simple way (equa- 
tion [5|, at least for barriers which are close to constant or 
square-root. Since we are in this regime, we assume that our 



prescription (i.e., C(m,t) = g(m)c(t\m)) works well enough 
for our purposes, even when the initial power spectrum is 
not white noise, and the barrier associated with the ran- 
dom walk problem which gives the mass function n that 
is of interest, is not constant. The second property is that 
this prescription reduces to the known exact creation rate 
in the white-noise case. This fact, while seemingly obvious, 
was never imposed as a require ment to be obeyed by t he 
creation rate in previous works (IBlain fc Longair|[l993a1 lb1: 
ISasakilll994l ; iKitavama fc Sutdll99d V 

Figure [7] compares this assumption with the measured 
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Figure 8. Conditional creation rates. Symbols, line styles and choice of masses as in Figure [6] 



creation rates in the simulations. The simulation measure- 
ments in Figures [7] and |4] are the same, except that in the 
latter, data is normalised in time. Notice that the heights 
of the curves increase with decreasing mass. This reflects 
that fact that, in hierarchical models, more small haloes are 
created during the history of the Universe than are massive 
halos (which are only created at later times). The constant 
barrier model works well for massive haloes, but it overpre- 
dicts the creation rate of less massive haloes - showing a 
similar discrepancy as in the Press-Schechter mass function 
in that mass regime. In all cases, the square-root barrier and 



the lSheth fc Tormenl (|l999l ^ creation rates match N-body re- 
sults reasonably well. 



4.2 Conditional rate 

We assume that the same prescription can be applied for the 
conditional case. I.e., the rate of creation of haloes of mass 
m at time t conditioned to belong to M haloes at a later 
time T > t is given by 

C(m,t \M,T) = pmN(m,t \M,T) J^J (17) 
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<|Shethll2003l ). Figure [5] shows our results for the same set of 
m and M as in Figure [S] As in the unconditional case, the 
simulation measurements in Figures [8] and [6] are the same, 
except that in the latter, data is normalised in time. 

For both choices of m, the height of the curves decreases 
with increasing final mass M (compare top and bottom pan- 
els). This effect is less pronounced for the smaller m (com- 
pare left and right panels). In all panels, curves based on 
the square-root barrier provide a more accurate description 
of the measurements. 



5 DISCUSSION AND CONCLUSIONS 

We have used the excursion set approach to study how the 
distribution of halo creation times is modified if haloes are 
assumed to form from an ellipsoidal rather than a spheri- 
cal collapse (see Figure [lj . The creation time distribution 
governs the time dependence in the creation term found 
in Smoluchowski-like interpretations of the e v olution of the 
halo mass function l|Sheth fc Pitman] Il997l : iBenson et al.l 
2005; Benson 2008) (i.e., normalising the creation rate in 
time gives the creation time distribution). In this approach, 
halo abundances and creation times can be derived from the 
study of random walks which cross barriers of some spec- 
ified height. The barrier itself is specified by the physics 
of gravitational collapse; that associated with spherical col- 
lapse has a constant height and this allows simple ana- 
l ytic solutions for halo abundances and crea tion times (e.g., 
feond et al.1 Il99ll : IPercival fc Miller! Il999h . Dropping the 
spherical collapse assumption is not trivial, as it results in a 
barrier whose height increas es with distance along the walk 
l)Sheth. Mo fc TormenluOOll ). Moveover, the moving barrier 
associated with ellipsoidal collapse is not well-suited to the 
study of creation times (see equation [2] and related discus- 
sion, as well as discussion in lSheth fc Tormenll2002j ). 

For this reason, we approximated the ellipsoidal col- 
lapse barrier with a square-root barrier (equation [<J]) for 
which an analy tic solution for the halo mass function exists 
|Breimadll966h . Moreover, this barrier yields a halo mass 
function which is in good agreement with simulations (Fig- 
ure [3]). We used Monte-Carlo realisations of random walks 
to show that the associated creation time distribution is re- 
lated to the halo mass function in just the same way that 
it is for the constant barrier (equation [5j) . Because the halo 
mass functions associated with the two barriers are differ- 
ent, the creation time distributions also differ (Figure fj}. 
The moving barrier based predictions were in slightly bet- 
ter agreement with measurements of halo creation times in 
simulations (Figures [4] and [5]) . 

We also presented the first derivation of the conditional 
creation time distribution using the moving-barrier model 
of ellipsoidal collapse. In this case, the differences between 
the constant and moving barrier predictions were somewhat 
more dramatic; the moving barrier predictions were in sub- 
stantially better agreement with the simulations (see Fig- 
ure [6] - especially the panels on the right) . For both the 
unconditional and conditional versions we showed there was 
a simple link between c and / (equations IS1 and 1 1 Of) . but 
showed that this link is special to the constant and square- 
root barriers. While not true in general (see Appendix A for 
a counter example), we argued that, for the mass functions 



of current interest in cosmology, this link should provide rea- 
sonably accurate approximations to the creation time distri- 
bution. 

We also presented an approximation for the normalisa- 
tion constant which converts the creation time distribution 
into a halo creation rate. This was motiv ated by connect- 
ing o ur results to coagulation theory (e.g. ISheth fc Pitman! 
1997). Figures [7] and [H] indicate that this yields reasonably 
accurate results. 

An alternative approach is to express the creation rate 
in terms of R(m, m'\t), the rate of mergers of m- haloes with 
m'-haloes, creating (m + m')-haloes as a result. This merger 
rate has been measured rec ently in cosmological simula- 
tion s dFakhouri fc Mal[2008aH bh and st udied in the spheri - 
cal (jLacev fc Colelll993h and ellipsoidal l)Zhang et al.ll2008al ) 
collapse versions of the excursion set approach. The prob- 
lem in this case is that, with the exception of white-noise 
initial conditions, the exc ursion set theory pre dicts that 
R(m,m'\t) ^ R(m\m \t) ijSheth fc Pitman! Il997h (but see 
iNeistein fc Dekell2008l . for a possible solution). Even if the 
merger rate were defined unambiguously within the excur- 
sion set approach with moving barriers, care must be taken 
when the integrals in equations (|12p and (|13p are com- 
puted. Our prescription of the creation rate circumvents 
these problems altogether. A third possible metho d is to use 
the merger ra te developed bv lBenson et al.l (|2005h (see also 
lBensorjf2008h . This method avoids the complications in the 
excursion set theory by estimating the coagulation kernel K 
numerically. Unfortuntely, the answer in this case depends 
on the choice of regularisation technique. The formulation a 
complete model for halo mergers with general initial condi- 
tions and with the advantages of ellipsoidal collapse remains 
an open, and quite interesting, problem. 

Halo creation rates have been used to model star 
for mation and AGN activi ty. We refer the reader 
Haehnelt fc Reed (Il993 



to 



Haim an fc Loebl dl99sh. iHaiman fc Menoul 
iHosokawal j|2002l) , [Cranato et alj (120041 ). lLapi et all 



Haehnelt ct al. 




and IWang et al.l 1 20091 ') for applications of halo creation 
where our results can provide an improvement. In some 
cases, the derivative of the mass function dn/dz is used 
to replace the creation rate - while in others, the creation 
rate is extracted from dn/dz without proper justification 
(see the references in the Introduction). It is our hope that 
our analysis will be useful for such studies. We caution, 
however, that a halo may be created smoothly through a 
process of gradual accretion, or violently, through a merger. 
If this difference matters, it may be more appr o priate 
to use a full merger history tree. iMoreno et. al. I (|2008h 
describ e s such moving-barr i er bas e d trees - in this context , 
also see Ifflotelis fc Popold (1200611: IParkinson et all l|2007h : 
INeistein fc Dekell i|2008l ) and lZhang et al.1 ^OOSbl 1 ). 

The phenomenon called 'down-sizing' is often associ- 
ated with two trends: (a) the tendency for the most massive 
galaxies to host the oldest stars, and (b) the tendency, at 
later times, for star formation to occur in haloes of lower 
mass. Our creation time distributions conditioned on final 
halo mass have been used to understand trend (a): if star 
formation only occurs in sufficiently ma ssive haloes , then 
(a) arises natural ly in hierarchical mod els ijSh eth 2003). This 
was confirmed bv lNeistein et all l)2006h . who then found that 
they were unable to explain trend (b) using this same mech- 
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anism. However, a corrollary of Sheth's argument is that, 
if star formation does not occur in haloes above a certain 
critical mass, then trend (b) is the result. That is to say, if 
stars only form in halos between a minimum and maximum 
mass range, then down-sizing trend (a) is the result of the 
minimum mass, and trend (b) is the result of the maximum 
mass. Therefore, we expect our results to provide further in- 
sight into this phenomenon in general, and into the critical 
mass scales in particular. 
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APPENDIX A: EXCURSION SET RESULTS 
Al Unconditional formulae 

In this section we provide expressions for f(v)dv and 

• Constant barrier (Press-Schechter) : 

f(u)du = x ffe-^. (Al) 
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• Square-root barrier (exact): 

w v 



2 /a(-/3) 



where D\(x) are the parabolic cylinder functions, 
= />(-/?) = f°°dxD 2 ( x ), 



(A2) 



(A3) 



and the {A}-eigenvalues satisfy D\(—[3) = (lBreiman|[l96r3 : 
iMahmood fc Raieshll2005l : iGiocoli etal]|2007M ~ 
• Square-root barrier (series approximation): 

f(u)Au= ,/^e" fll+ ^ ]2 [l + ^] — , q~ 0.2461. (A4) 

V 27T a/!/ 1/ 



See lSheth fc Tormenl (|2002l ). 
• 'Sheth-Tormen' result: 

f(v)dv = A x [^-e-^[l + v- p ] — , (A,p) = (0.332,0.3). 

y 2lT V 

(A5) 

See lSheth fc Tormenl (Il999l ). 

In the above models we have suppressed the parameter q 
of equation @, which is (1, 0.55, 0.55, 0.707) respectively for 
each model. This can be incorporated by replacing v — > qv 
in the corresponding formulas. 



A2 Creation times from Bayes' rule 

Consider the joint probability that a random walk first up- 
crosses a constant barrier (with 5-intercept between S c and 
8 C +d<5 c ) between S and S' + dS'. Using Bayes' Theorem, this 
can be written as 

P(S,S c )dSdS c = f(S\5 c )f(5 c )dSd5 c = c(S c \S)c(S)dS c dS, 

For constant barriers, IPercival fc MilleU l|l999l ) argue that 
the 5 c -prior must be uniform. First they note that all walks 
must have a creation event for any barrier, regardless of 
its height (see Figure [T] left panel). Moreover, for any two 
equal-sized intervals dS c i and d8 C 2, the probability that such 
a creation event exists must be equal. This is because the 
steps in the walk are uncorrelated, implying that any point 
along the walk can be regarded as the starting point of a new 
walk. Therefore, the walk is not altered at different values 
of 8c and the probability of crossing two different barriers at 
some point is the same. In other words, the <5 c -prior is given 
by 



(A6) 



where the constant A8 C is infinite, since S c £ [0, oo). 

Given J(S c )d5 c and f(S\8 c )dS, one can marginalise the 
joint distribution in S c . This yields 



c(S)dS = dS / d<5 c [P(S,<5 c )] = 



dS 



2AVSAS C 



(A7) 



where A = (y^Tv/2, 2, 2.08, 1.893) respectively for each of 
the four models presented above. Inserting (|A6|I and ()A7|) 
in Bayes' formula gives 



c(5 c |S)d5 c = 2AVSf{S\Sc)d5 c 



But f(S\6. 
(2u c /S c )c(u, 



= {y c /S)f{y c 

C ), from which 



{S is fixed) and c(<5 c |S) 



c(^c)d^ c = Ayfv~cf(yc)dvc 



(A9) 



Thus, for constant barriers, there is a remarkably simple 
relation between the creation distribution c and the first 

cros sing distributio n f. 

IPercival et all |2000) argue that equations (IA6|I 
and (|A9[) remain true for all choices of /, which in the 
present context mean, for all barrier shapes. We will now 
show why this is incorrect, and also why equation ()A9[) nev- 
ertheless provides a good approximation to the correct an- 
swer. 



A3 Why it doesn't work in general 



Consider the linear barrier 



B{S,5c)=8 c [l 



(A10) 



This is special case of the barrier in equation |2]) with 7 = 1 
(we have suppressed the parameter q). Notice that we have 
deliberately replaced f3 — ► — /3 in the above expression. In 
this work we will only consider the j3 > case to avoid 
barriers that intersect. Moreover, if j3 > 0, all walks are 
guaranteed to cross. The exact solution is known, and it is 
given by 



f(u)du 



.| ( i_£)3di/ 
v 



(All) 



where v = 81/ S (fixed S) (|Schroedingerlll915l : IShetblll998l ). 
If equation (|A6(I holds, then the creation time distribution 
should be given by 



c(vc)dv c = A^fv~cf(yc)dvc, where A ■- 



TV _ 

2 6 



(A12) 



(A8) 



and v c — 8c /S (fixed 8 C )- 

To test this, we performed a Monte Carlo simulation of 
the mass histories associated with random walks where lin- 
ear barriers (with f3 = 1) are used to select creation events. 
Figure I All shows one sample mass history of this ensemble. 
Compare this with Figure [1] which shows the same process, 
but with constant and square-root barriers. The jagged line 
is a random walk, and the brown solid circles are the as- 
sociated history. The long dashed lines denote jumps in the 
history. For completeness we have included the linear barrier 
with 8 C — 8co, depicted as a dotted line in the lower left. Fig- 
ure E2] shows our Monte Carlo data (brown open squares). 
The solid black curve is the prediction in equation (|A12|) . 
This disagreement invalidates the claim that c(i>c)dvc is al- 
ways proportional to •Jvc'j{v^)dvc- 

Before moving on, note that the Monte Carlo data in 
Figure IA12I follow a smooth curve that is quite different 
from that associated with constant or square-root barriers 
(Figure E| . The main difference is that the latter peak at 
some intermediate value of v c , whereas the former decreases 
monotonically with v c - This is because it is unlikely for a 
random walk to upcross a constant barrier (or a square-root 
barrier) after a few steps. In other words, creation events 
with S <C 8c (i.e., v c ^ 1) are very unlikely. Similarly, it 
is unlikely that a walk survives for many steps without be- 
ing absorbed by a constant barrier. That is, creation events 
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5 10 15 20 

S 

Figure Al. The mass history associated with a random walk 
(jagged line) and linear barriers. The brown filled circles on the 
random walk denote the creation events. Notice that the barriers 
become steeper as S c increases, and the separation between any 
two barriers increases with increasing S. For reference, the dotted 
line in the lower left represents the barrier associated with <5 C = 
<5 c o- 



i mm 1 — i — i i i i 1 1 1 1 — i — i i i i 1 1 

Bayesian result 




0.1 1 10 100 

Figure A2. The creation time distribution associated with lin- 
ear barriers in self-similar form. The variable is c denotes at 
fixed mass. The squares show the distribution measured from an 
ensemble of random walks with linear barriers. The solid curves 
show the associated predictions - assuming that equation l|A6|l 
applies (the 'Bayesian' result). The discrepancy between the the- 
ory prediction and the data indicates that care must be taken 
when using the result in equation $5\ . 



with S 3> 8% (i.e., v c 3> 1) are unlikely. For the linear bar- 
rier, the ensemble of creation events is dominated by points 
with small S c (i.e., small This is because the linear bar- 
rier becomes steeper as S c decreases. Since the height of a 
linear barrier decreases with S, it becomes easier for walks 
to upcross a barrier as S c decreases. This even allows for 
cases where creation events are selected from points along 
a walk with 8 < 0. Such cases would be impossible for the 
constant and the square-root barrier models. 

In principle, the creation time distribution with lin- 
ear barriers can be computed analyticall y, without using 
the Bayesian approach presented here (e.g.. lKarlin fc Tavlorl 
Il975l ). Since the mass function associated with linear bar- 
riers does not resemble halo abundances in N- body simula- 
tions, we do not pursue this any further (but see lShethlll998l . 
for interesting applications of this barrier). 

A4 Why it is a useful approximation in practice 

If equation (|A9(I is incorrect in general, then why then did it 
work so well in the main text? The first step is to recognize 
that, because of the property highlighted by equation @, 
a uniform distribution in S c is appropriate for the family of 
square-root barriers of interest to us (equation [3]) , and so, 
for square-root barriers, equation (| A9f) is exact. Our Monte 
Carlo simulations shown in the main text confirm that this 
is indeed the case. 

The second step is to note that, for barriers with gen- 
eral 7, the difference between two barriers carries additional 
factors of 5 C (e.g., equation[2]). As a result, assumption (IA6|I 
is no longer valid. E.g., for linear barriers, 

B(S,S c2 )-B(S,S cl ) =6 c2 -6 cl -p(^L -■!?-) ^802- Soi. 

V0 C 2 del/ 

(A13) 

This property makes the linear barrier considerably different 
from the constant and square-root barriers. Figure |X2] shows 
that the distance between any two linear barriers increases 
with increasing S (compare with Figure [TJ, so the <5 c -prior 
is not uniform. 

It is important to emphasize that the central conclusion 
of this and the previous subsection - that the assumptions 
behind equation (JSJ) do not hold in general - do not depend 
on the fact that the linear barrier decreases in height with 5*. 
E.g., barriers of the form given by equation @ with /3 > 
and < 7 < 1/2 increase with S. However, 

B(S,8 c2 )~B(S,8 cl )^S C 2-8 cl ; (A14) 

the separation between any two barriers increases with S, 
so the (5 c -prior is not uniform in this case either. For these 
barriers too, equation (1 A9f) is incorrect. Nevertheless, for 
7 close to or 1/2, equation (|A9[) should provide a rea- 
sonable approximation. This is the fundamental reason why 
barriers with 7 = 0.6, or of the form required to give the 
ISheth fc Tormenl l|l999h formula as the first crossing distri- 
bution, are likely to have creation time distributions which 
are well approximated by equation (|A9|I . 

A5 Conditional formulae 

The conditional mass function and creation times distri- 
bution can be obtained from the unconditional crossing 
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distribution of the barrier B in equation JB). The natural 
variables in this case are s/So and 77,3 = a c /vSo (where 
a c = Si — So — f3\fSo). Explicity, the crossing distribution 
can be written as 



f(s/So\vp)d(s/S )=g(s/S ,r,p) 



d(s/S ) 
s/So + 1 ' 



(A15) 



where two forms of g are given below (|Breimanl 1 19661 ; 
ISheth fc Tormenll2002h . 



• Square-root barrier (exact): 



g(s/S ,r)f}) = 



{A} 



(s/S + 1) A / 2 ' 



(A16) 



• Square-root barrier (series approximation): 



g(s/So,V/3) 



\Vp + /VVS> + l[l + a (s/So) ] J 
s/2ns/S 

r (vp + PVWo + i) 2 1 s/So + 1 

XeXP l 2s7So 1"" 



(A17) 



where 



*(s/S ) 



^s/s?TT' = and "» = (i-^K-i- 



E 



• B ayes' rule: 

The joint distribution of s/So and 77/3 is given by 

P(s/S ,r ] p)A(s/So)dr ! p = f(s/So\vp)f(vp)d(8/So)dvp 

= c(r ! p\s/So)c(s/So)d(s/So)dr 1 p.(A18) 



Following equation (|A6|) for Si and So and using the fact that 
So is fixed, it can be shown that the 77/3-prior is uniform: 



(A19) 



where A77 is an infinite constant. Marginalising over the joint 
distribution in r\p we obtain 



c(s/S )d(s/So) 



G d(s/S ) 
s/So + 1 At? ' 



where 



G(s/S )= dr)p g(s/ So, rip). 



(A20) 



(A21) 



APPENDIX B: COAGULATION THEORY 
Bl White- noise initial conditions 

The creation rate in the discrete Smoluchowski equation is 

C(m,t)= 22 — — n(m \t)n(m-m \t). (Bl) 

— 1 

If the kernel A' is additive, the solution approaches the 
white-noise Press-Schechter mass function in the continuum 
limit (large m and small S c ). Moreover, ISheth fc Pitman! 
(1997) showed that the above equation can be written as 



C(m, t) = n n(m\t) 



m — 1 1 d5 c 



l + S c \ At 



Pr 



m ' \ m — m' • 



(B2) 



where n is the mean number-density of particles. If we 
picture haloes as a collection of particles held together 
by (m — 1) non- intersecting bonds (branched polymers), 
Pm'lm-m' gives the probability of obtaining an m'-halo and 
and (m — m')-halo by deleting one random bond in the m- 
halo. The sum over a normalised probability is trivial, so the 
creation term is simply 

„, > _ , . . m — 1 1 dS c 1 
C(m,t) — nn(m\t)- 



1 + 5, 

In the continuum limit, this becomes 
C(m,t) = pmn(m\t) 
Similarly, in the conditional case, 
C(m,t\M,T) = nN(m\t,M,T) 



At 



d4 
At 



1 + 5, 



l\d5 c 
dt 



l|Shethll2003h . In the continuum limit, this becomes 
C(m,t\M,T) = pmN(m\t,M,T) |^|. 



(B3) 



(B4) 



(B5) 



(B6) 



Notice that the two rates (equations IB4I and IB6|) are 
related by the following consistency condition: 



f 



AM C(m, t \M, T) n(M\T) = C(m, t). 



(B7) 



In other words, the unconditional rate is recovered from the 
conditional rate by multiplying by the number density of 
M-haloes at T and integrating over all possible M > m. 
This consistency relation is true in general, not just when 
the initial conditions are white-noise. 



Inserting equations (|A15|) . (|A19I) and (|A20|) in Bayes' rule, 
we obtain 



c(r)g\s/S )dr) S = ^dr/p. 



Comparing (|A15() to (|A22|) . we find that 

c(rip\s/S )drip = ^ (s/So) /(s/So|^)dr/ /3 , 

where 

As/ s 0) = (s/So + 1)/G. 



(A22) 



(A23) 



(A24) 



This result is confirmed by our Monte Carlo simulations 
(Figure El). 



B2 The time-normalised creation rate 

In this appendix, we show that creation rate C(m,t) (equa- 
tion !14l Section 4) is related to the creation time distribution 
c(t\m) (equation O Section 3) in a simple way. First, notice 
that c(t\m) can be written in terms of v c = 5^/S. That is, 



c(t|m)=c(5 c |S)|^| 



. , I d^ c 1 1 AS C I 2i^ c I A5 C I 



(B8) 

where we have used the fact that |d In v c /A In <5 C | = 2. 

Now, following equation (0, the above expression can 
be written as 



c(t\m) = y^^AA^/C^c)^! 



(B9) 
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On the other hand, 

where we have used the fact that |dln v/dlnS\ = 1. 

Our prescription for the creation rate (equation 1 14p is 

_ 2 mn(m\t)\dS c \ _ 2 v | dS 1 1 d<5 c I 
Taking the ratio of (|B11[) and (|B9(I . we obtain 



C(m,t) 


f 


s/w|n| 


1 dS c 1 
dt 


P 2 1 d5 1 


c(t\m) 


5c 




1 dt 1 


2^\/5 dm 



(B12) 

where we have dropped the distinction in notation between 
v and v c . Simply put, 

C{m,t) = g(m)c(t\m). (B13) 

Integrating both sides with respect to time, 

j C(m,t)dt = J g{m)c{t\m)dt = g(m) J c{t\m)dt = g(m), 

(B14) 

since c(t\m) is, by definition, a time-normalised distribution. 
Thus, 

, . , C(m,t) C(m,t) ,„,_v 
c(tm) = \ '.' = , y (B15) 

In other words, c(t\m) can be obtained from C(m,t) by nor- 
malising the latter in time. 
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